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ABSTRACT 

' We present a first cosmological analysis of a refined cluster catalog from the Red- 

Sequence Cluster Survey (RCS). The input cluster sample is derived from 72.07 square 
degrees of imaging data, and represents a deep uniform subset of the imaging data 
available within the RCS survey images which probes to the highest redshift and low- 
est mass limits. The catalog contains 956 clusters over 0.35 < z < 0.95, limited by 
cluster richness and richness error. The calibration of the survey images has been ex- 
tensively cross-checked against publicly available Sloan Digital Sky Survey imaging, and 
the cluster redshifts and richnesses that result from this well-calibrated subset of data 
are robust. We analyze the cluster sample via a general self-calibration technique in- 
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OO . eluding scatter in the mass-richness relation, and using reasonable priors on the Hubble 

. constant. We fit simultaneously for the matter density, the normalization of the 

O ', power spectrum, as, and four parameters describing the calibration of cluster richness to 

I mass, its evolution with redshift, and scatter in the richness-mass relation. The princi- 

I pal goal of this general analysis is to establish the consistency (or lack thereof) between 

Qh' the fitted parameters (both cosmological and cluster mass observables) and available 

Q ■ results on both from independent measures. From an unconstrained analysis, the de- 

-4— » ' rived values of Jim and erg are 0.31]'^q']^J and 0.67l[J;{^ respectively. An analysis including 

. Gaussian priors on the slope and zeropoint of the mass-richness relation gives very sim- 
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ilar results: 0.30_q']^^ and 0.70_g'^5. Both analyses are in acceptable agreement with 
^ I the current literature. The derived parameters describing the mass-richness relation in 

H ' the unconstrained fit are also eminently reasonable and in good agreement with exist- 

ing follow-up data on both the RCS-1 and other cluster samples. Our results directly 
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demonstrate that future surveys (optical and otherwise), with much larger samples of 
clusters, can give constraints competitive with other probes of cosmology. 

Subject headings: galaxies: clusters: general, cosmological parameters, methods: data 
analysis 

1. Introduction 

The dependence of the number density of massive clusters with redshift on the cosmological 
mass density, O^, and the normalization of the power spectrum, erg, has been noted in the literature 
for nearly two decades (e.g., Evrard 1989; Oukbir & Blanchard 1992; Fan, Bahcall & Cen 1997) 
and cosmological tests have been attempted using cluster counts over a similar span (e.g., Edge et 
al. 1990; Oukbir, Bartlctt & Blanchard 1997; Eke et al. 1998; Henry 2000; Bahcall & Bode 2003; 
Pierpaoli et al. 2003; Henry 2004). Convergence on the values of O^, and particularly on ag has 
been slow using this general method, though some recent results tend to favor a high normalization, 
low-density model (e.g., Bahcall &; Bode 2003), irrespective of cosmological constraints from other 
techniques (e.g., see Spergel et al. 2003, and references therein). 

Recently, various extensions of these techniques have been suggested (Haiman, Mohr &; Holder 
2001; Levine, Schulz & White 2002; Hu & Kravstov 2003; Hu 2003; Majumdar & Mohr 2003) as a 
method for measuring the equation of state parameter of dark energy, w, and a number of authors 
have presented parameter accuracy estimates based on future surveys (e.g., Kneissl et al. 2001; 
Weller et al. 2002; Newman et al. 2002; Weller & Battye 2003; Majumdar & Mohr 2004; Wang et 
al. 2004) . Despite this flurry of recent work, little progress has been made in actually confronting 
the suggested techniques with real cluster survey data, and it remains to be seen whether systematic 
issues or unconsidered physical effects will limit the utility of cluster mass function cosmological 
measurements. 

In this paper we present a first analysis of a completed cluster survey - the Red-Sequence 
Cluster Survey (RCS-1; Gladders & Yee 2005, hereafter GY05) - using the self-calibration technique 
suggested by Majumdar & Mohr (2003). Unlike previous analyses of X-ray data (e.g., Bahcall & 
Bode 2003; Pierpaoli et al. 2003), the cluster sample used here is drawn from a single homogeneous 
survey which was executed precisely to enable this type of cosmological analysis, and contains a 
much larger number of clusters spanning a broad range of mass and redshift. 

This paper is organized as follows. In §2 wc describe the input cluster catalog, and exten- 
sive tests designed to demonstrate that it is robust and well calibrated. Section 3 describes the 
methodology of our self-calibration analysis of this catalog, which simultaneously constrains both 

the cosmology and the mass-observable relation used. In §4 we describe and discuss the results 
in the context of current cosmological parameter constraints and our understanding of the mass 
richness relation in clusters. Our conclusions are summarized in S5. 
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2. Cluster Catalogs 

The entire RCS-1 survey contains a total of ~ 90 degrees of imaging data from both the Canada- 
Prance-Hawaii Telescope (CFHT) and the Cerro Tololo Inter-American Observatory (CTIO) 4m 
telescope with coverage in both the Rc- and z'-band. Exhaustive details of the data acquisition, 
reduction and analysis and the generation of cluster catalogs from these data are presented in 
GY05. Some of the data were taken under sub-optimal conditions, with significant transparency 
and seeing variations, compromising both the image depth and the photometric stability. Based 
on records from the observing runs, and an examination of the reduced images and catalogs, we 
have eliminated the poorest data, and restrict the analysis which follows to a survey area of 72.07 
square degrees. 

The primary parameters of each cluster relevant for the cosmological measurement performed 
here are the estimated photometric redshift and the cluster richness, which both rely on stable, well- 
calibrated photometry for accuracy. In order to check the RCS photometry we have exhaustively 
compared the portion of the RCS photometric database which overlaps with the most recent Sloan 
Digital Sky Survey (SDSS; Adelman-McCarthy (2005)) public data release. The overlap consists 
of 26.29 square degrees, of which 1.44 square degrees is in the CTIO portion of the RCS. We have 
applied extinction corrections derived from the atlas of Schlegel, Finkbeiner &: Davis (1998) to 
both the RCS-1 and SDSS magnitudes, and then checked the photometry for stability both at the 
camera level (i.e., inter-chip) and from pointing-to-pointing and run-to-run. Comparisons in the 
z'-band have been made directly; in the i^c'-band we adjust the zeropoints to eliminate the scatter 
in comparison to the r'-band, while preserving the mean offset. We have checked and corrected the 
RCS-1 color terms for each mosaic camera, as well as the zero points for each chip. The corrections 
are typically a few hundredths of a magnitude, the expected value given the nominal photometric 
calibration of the survey (GY05) . 

We have used the bright star and bright galaxy color distributions from each pointing in the 
RCS-1 to check for color offsets between pointings, in a method similar to that of Hsieh et al. (2005). 

The comparison of our derived color offsets to direct color corrections derived from comparison to 
SDSS magnitTidcs shows this provides accurate colors to better than 0.03 magnitudes. Internally, 
the color corrections for each pointing deduced from the star and galaxies separately have a scatter 
of less than 0.025 magnitudes. 

The color corrections described above are sufficient to ensure accurate redshifts for RCS- 
1 clusters. However, the richness measurements require well measured magnitudes, and so we 
have developed a relatively simple correction scheme which uses the counts of faint galaxies in 
combination with the color corrections. To the precision we require, galaxy counts per pointing 
should be constant, both because of the large area of each pointing, and because the imaging is 
sufficiently deep to probe a large redshift column. A comparison of this simple approach in the 
regions of SDSS overlap shows that it approximately halves the residual magnitude errors, to about 
0.027 magnitudes in each filter per pointing; we have applied this method to all of the data discussed 
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here. Further details of this method will be provided elsewhere. 

After the procedure described above, the final photometry for the best 72.07 square degrees 
of the RCS-1 is well calibrated to within a few hundredths of a magnitude. This is sufficient 
to provide accurate redshifts, with typical uncertainties of less than 0.05 at z < 0.6, increasing 
to about 0.09 at z = 1 (e.g., Gladders 2003), and accurate richnesses (Yee & Lopez-Cruz 1999). 
The resulting stability of the magnitudes and counts of galaxies across the entire dataset ensure 
that the richnesses are systematically stable to well below the random uncertainties which arise 
predominantly from shot noise in the galaxy counts. The final cluster catalogs based on these 
data thus have both accurate redshifts and richnesses - precisely the data required to make a 
cosmological measurement. 



3. Self-Calibration Analysis 

The input to the self-calibration analysis is a catalog of RCS clusters down to a 3.3-sigma 
significance (GY05) limit. As detailed in GY05 and Gladders & Yee (2000), this significance is 
determined from a detailed bootstrap analysis of the non-cluster regions of the RCS-1 survey data 
directly. At the our chosen 3.3-sigma limit, the contamination of the sample is less than 5% based 
on our current understanding of the RCS-1 false-positive rate (Gladders 2001). As in GY05 the 
nominal lower redshift limit of the catalog is z = 0.2. However, to avoid a possible Malmquist bias 
from clusters scattering out of the sample to below the redshift limit of the catalog, we limit our 
analysis here to clusters at z > 0.35 (i.e., more than 3-sigma above the lower redshift limit of the 
catalog). We also explicitly limit the catalog to z < 0.95; the detection limit changes rapidly at 
z > 1 as the 4000 A break moves into the z'-band, and clusters at the highest redshifts require more 
careful calibration of the photo-z and richness measurements before they can be reliably used for 
cosmological tests. 

We include in our analysis our current understanding of the RCS selection functions (Gladders 
2001) with respect to both richness and blue fraction. We take the evolution of the RCS-1 cluster 
blue fraction with redshift from Loh ct al. (2006, in preparation). For all redshifts considered 
here the completeness corrections arc less than 25%, and significantly less than that in all but the 
highest redshift bin. We have not attempted to model the error on the derived selection functions, 
leaving that effort to a future paper. Given the magnitude of the corrections however, we do not 
expect uncertainty in the catalog completeness to be a dominant source of error 

We use the richness parameter Bgc (Longair &; Seldner 1979) as our estimator for cluster 
richness. Bgc is the amplitude of the cluster-center to galaxy correlation in units of {h^^ Mpc)^'^^ 
(the fundamental measurement is the evolution-corrected, background-corrected and luminosity- 
function normalized overdensity of galaxies within a fixed proper radius of 500 kpc, with h = 0.5). 
Bgc has been shown previously to correlate with cluster mass (Yee & Ellingson 2003; Yee & Lopez- 
Cruz 1999). We specifically use the red-sequence richness, BgcR (GY05), in all analyses which 
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follow; BgcR is expected to more closely trace the evolution of mass since red cluster galaxies are 
well established even at z = 1, and the measurement errors on BgcR are much smaller than on the 
total richness. In addition to the significance threshold, we limit the input catalog to only clusters 
with BgcR > 300, and further limit the catalog to clusters where the error on BgcR is less than 50%. 
The BgcR > 300 cut ensures a monotonic mass limit with redshift; the significance cut is insufficient 
to achieve this since the RCS-1 data are more sensitive to clusters at moderate redshifts. The limit 
on richness errors serves to eliminate clusters which are poorly measured due to having only small 
numbers of galaxies (usually because they are near some artificially shallow portion of the survey 
data - such as near a bright star). In total the resulting catalog, including both the redshift and 
richness limits described above, contains 956 clusters. The total number of clusters represented by 
the catalog, once corrected for incompleteness, is 1086. 

The expected surface density of clusters in a solid angle AO at redshift z to a limiting mass 
Miijn from a fiducial mass function ^ is 

CO 

— {z) = AO— — (z / -rrjdM . (1 
dz ^ ' dzdn^ ' J dM ^ ' 

Assuming the data are grouped into redshift bins of width Az, the directly observable quantity 

z+Az/2 

N{z) will then be given by N{z) = J ^{z')dz'. The adopted form of the mass observable 

z-Az/2 

relation (used to calculate Mn^), following Yee &; Ellingson (2003), is given by 

M200 = 10^^-i?,",«(l + z)^ (2) 

where 7 allows for any possible unknown evolution of the mass-richness relation. We use the Jenkins 
et al. (2001) mass function in our theoretical calculations and a simple NFW profile (Navarro, 
Frenk &; White 1997) for the dark matter halos to convert from M200 to MjcnHns- Any theoretical 
uncertainties in the conversion between masses are subsumed in the cluster mass-richness parameter 
constraints in our self-calibration analysis, described below. 

To compare theoretical predictions of the cluster redshift distribution to the RCS-1 input 
catalog and subsequently estimate the cosmological parameters, we use a Markov-Chain Monte- 
Carlo (MCMC) analysis (S. Majumdar & G. Cox 2006, in preparation). It has already been 
demonstrated (Levine, Schulz & White 2002; Majumdar k Mohr 2003, 2004; Lima k Hu 2005) 
that cluster surveys with thousands of clusters would have enough information (in theory) to 
simultaneously determine both cosmology and cluster physics, providing a direct measure of the 
cluster mass-observable relation in addition to cosmological constraints. 

Our basic set of parameters consists of 3 cosmological parameters (Omj^s, and h) and 4 
cluster parameters. These are the amplitude Asgc and slope a of the mass-BgcR relation, its 
redshift evolution parameter 7, and the fractional scatter fsc in the m.ass-BgcR relation. Due to 
the exponential sensitivity of the cluster redshift distribution to the underlying cosmology, the 
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cluster counts are not only sensitive to the mean richness-mass relation (given by Equation 2) 

but also to the actual distribution including scatter. We have modeled this by a gaussian scatter 
parametrized by the mass independent fixed scatter fraction fsc- Wc do not have any implicit 
redshift dependence on the scatter. The scatter is incorporated by multiplying the mass function 
with a "selection function" F{M, z) such that 

With the inclusion of the selection function, the lower limit of the integral over mass is changed 
from Miirn{z) to some Miow{z)- For our calculations, we have fixed M;o^ = 8 x lO^^M©. For 
our best fit M/j^, the final constraints are not too sensitive to the value of Mi^y, as long as it is 
< IO^^Mq. However, we have seen that with the inclusion of scatter one needs to increase the 
accuracy of numerical integral routines to get convergent 

Note that from the consideration of cluster counts, the amplitude Asgc and slope a are de- 
generate (by construction) and the cluster redshift distribution only really constrains the limiting 
mass Miimo and its redshift evolution. Hence, we present the basic analysis which follows in terms 
of Mii„,^y Results involving Abqc and a arc given for comparison with direct cluster observations. 
This degeneracy between Abqc — a is broken when one takes priors (or a joint analysis) of targeted 
observations which provide calibrating constraints on cluster mass. A more complex analysis than 
that presented here which fits the mass function with redshift (rather than its integral) would also 
in principle break this particular degeneracy. However, such an analysis is not possible with the 
limited size of the RCS-1 survey. 

In our analysis we assume a flat ACDM Universe (i.e., w = —1). Cluster surveys alone 
cannot constrain the Hubble constant and so we put a Gaussian prior on the Hubble constant (i.e., 
h = 0.72 lb 0.08). We also fix the spectral index Ug and the baryon density to first year WMAP 
values (Spergel et al. 2003). Other parameters are constrained by weak uniform priors, in order 
to constrain the parameter spaced searched in the MCMC analysis. These various priors are listed 
in Table 1. 

Our MCMC analysis uses the algorithm proposed by Metropolis et al. (1953) to randomly 
sample the parameter space with a Markov chain whose distribution asymptotically approaches the 
distribution from which it is being sampled (e.g., Lewis &; Bridle 2002). To construct the chain, 
we calculate the likelihood at each point in the parameter space under the assumption that the 
distribution of the clusters in a redshift bin is essentially Poissonian in nature. Our choice of redshift- 
bin thickness is optimal since the covariance along the z-direction is negligible for Az ~ 0.1 (Hu 
&; Cohn 2006). Moreover, the thickness is greater than the redshift uncertainties. Smaller redshift 
binning would require use of more generalized likelihood functions (Holder 2006; Hu &; Cohn 2006). 

Typically we run ~ 4 — 6 chains and need more than half a million points to reach convergence. 
We have checked that the chains span a large parameter space and sufficiently overlap with each 
other. We have also seen that inclusion of cluster parameters requires more time for the chains to 
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converge than those having no cluster parameters, since cluster parameters significantly widen the 
parameter space. 

4. Results and Discussion 

4.1. Analyses Without Mass- Richness Priors 

The principal goal of this paper is to provide a first observational test of whether large cluster 
surveys - in particular large optical cluster surveys - in combination with a self-calibration analysis 

can yield useful cosmological constraints. As such we arc interested in analyses with only weak priors 
on the fitted parameters (i.e., the utilitarian uniform priors listed in Table 1), as the consistency 
of the results to the literature (both cosmologically, and in the cluster-mass observables) yields 
important insight into the validity of this approach. However, we are ultimately interested in the 
best possible cosmological constraints from the existing RCS-1 survey data; and so we also include 
in the next subsection a preliminary analysis of the cluster sample using Gaussian priors on Abqc 
and a, derived from extant analysis of follow-up studies of the RCS-1 (Blindert et al. 2006, in 
preparation) and the CNOC-1 survey (Yee & Ellingson 2003). 

The results of our unconstrained analysis are shown in Figures 1-4. Figure 1 shows the input 
redshift distribution, and the best fit cosmological model. The MCMC analysis has been done 
in this case with M^jmo as one of the cluster parameters. Figure 2 shows the fully marginalized 
likelihood distributions for the six fitted parameters in this case, as well as the central values. These 
values, with 68% confidence limits, are summarized in Table 2. We show in Figures 3 and 4 the joint 
likelihood distributions for most parameter pairs. Figure 4 specifically focuses on the relationship 
between the fractional scatter, f^c, and other parameters. The most significant degeneracy, in terms 
of its impact on the cosmological results, is the relationship between the scatter in the mass-richness 
relation, and as- 

The final values of the cosmological parameters summarized in Table 2 agree well with recent 
results from the literature. In particular, our result on erg is in good agreement with the recent 
year-three WMAP constraints (ag values range from 0.722 to 0.772 depending on which datasets 
are analyzed in combination with WMAP; Spergel et al. 2006). Our value of similarly agrees 
with year-three WMAP (Spergel et al. 2006, ranges from 0.238-0.266), and spans results from 
combined analyses of WMAP and the SDSS (0.30±0.04; Tegmark 2004) and WMAP and the 2dF 
(0.231±0.021; Cole et al. 2005). 

Note that some recent cluster-based results tend to favor higher values of erg , which our analysis 
does not support. These results include studies of aggregate cluster samples selected by other means 
(e.g., Bahcall & Bode 2003) as well as the cluster-centric interpretation of the excess small scale 
power in the CMB in various experiments (Komatsu &: Seljak 2002; Goldstein et al. 2003; Bond et 
al. 2005). It has been already pointed out by a number of authors (Pierpaoli et al. 2003; Rasia et 
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al. 2005) that to get a better handle on as it is necessary to take into account the uncertainty in the 
scaling relations, and scatter. As pointed out by Seljak (2002), choice of the normalization of the 
scaling relation (either observational or from simulations) can give vastly different erg; an a priori 
choice of fixed scaling relation can give tighter yet biased constraints on cosmological parameters. 
However, a full marginalization over cluster variables is often not done. A comparison of our results 
to analyses using full marginalization of the cluster scaling relation shows excellent agreement. For 
example, Pierpaoh et al. (2003) ag = O.lltoM and Henry (2004) find erg = 0.66 ± 0.16 after 
marginalizing over the amplitude of mass-observable relation. In comparison to previous works, 
we not only marginalize over the normalization of the cluster scaling relation but also take into 
account any uncertainties arising from our incomplete knowledge of the slope, redshift dependence 
and scatter in the mass-observable relation by marginalizing over these parameters as well. 

It is also worth noting two additional points. First, the O^-ug degeneracies are weaker here 
by comparison to some methods and complementary to the stronger degeneracies seen in current 
X-ray cosmological constraints using clusters over a smaller redshift baseline (e.g., Pierpaoli et 
al. 2003), the CMB, or weak lensing (e.g., compare the ag-^lm panel of Figure 3 here to Figure 5 
in Van Waerbeke, Mellier & Hoekstra 2005). This is due to the long redshift baseline of the RCS-1 
cluster sample. Second, unlike a common tendency in the literature, we report full error bars on 
parameters (rather than, say, ug at a fixed value of 0^^=0.3). Encouragingly, the uncertainty in 
our constraints on cosmological parameters from the RCS-1 clusters are in excellent agreement with 
simple 'Fisher-Matrix' forecasts for upcoming large cluster surveys surveys. For, example typical 
surveys having 10000 — 20000 clusters arc predicted to constrain cg and to 0.07 and 0.03 
from dN/dz alone (Majumdar & Mohr 2004, Wang et al. 2004). RCS-1 with roughly 10 times 
less clusters gives errors on cosmology which are a factor of 3 larger. This consistency both in 
central values and uncertainties in cjg and Q.m implies that cluster surveys are on the right track 
and capable of giving us interesting constraints on cosmological parameters once more clusters are 
added to the analysis. 

The self-calibration approach also provides measures of the mass-observable relation and the 
evolution of its zero-point with redshift (see Equation 2), as well as the fractional scatter. To make 

the comparison to current data on the chistcr mass-richness relation, we have repeated the analysis 
of the RCS-1 cluster catalog by replacing M^jm,, with the degenerate parameters Asgc and a. We 

-1—0 07 _l_n Q1 .-M-1 

find values of Abqc = 10.55_]^'y]^ and a = 1.64_Q'gQ. These should be compared to the observed 
value for the CNOC-1 analysis X-ray selected clusters (Yee k, Ellingson 2003), which derived best 
fits values of Abqc = 9.89 ± 0.89 and a = 1.64 it 0.28 (corrected to h = 0.72, but not corrected 
for possible evolution as indicated by 7), as well as a recent analysis (Blindert et al. 2006, in 
preparation) of similar spectroscopy of 33 RCS-1 clusters at median redshift of z = 0.33, which 
yields similar values. Furthermore, the fractional scatter in mass observed in the Blindert et al. 
work is about 70%, in good agreement with our value of 0.73 it 0.22. 

That a generalized self-calibration analysis of the RCS-1 optical cluster catalog yields a mass- 
richness relation which is in agreement (in amplitude and slope) with a dynamical study of a 
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completely separate X-ray selected cluster sample (Yee & Ellingson 2003), and (in amplitude, 
slope, and scatter) with a dynamical study of an intermediate redshift subset of the RCS-1 clusters 
(Blindert ct al. 2006, in preparation), is a significant endorsement of the robustness and reliability 
of the RCS-1 optical cluster catalog, the self-calibration methodology, and the values of flm and 
(Tg derived from these cluster data. Moreover, we note that the mass-richness relation is in excellent 
agreement with a preliminary analysis of the weak lensing shear from RCS-1 clusters measured with 
the same survey data (Hoekstra et al., in preparation). 

Interpretation of the evolution term in the mass-richness relation is less obvious. Both data 
(e.g., see Lin &; Mohr 2004, and the discussion therein) and simulations (Zentner et al. 2005) 
suggest that as much as 50% of the light in cluster galaxies may be incorporated into the intra- 
cluster medium by z = (and hence would not counted in a richness measure if this intra-cluster 
light results from the destruction of significant numbers of bright galaxies). Additionally, an analysis 
of the X-ray temperatures of a small subset of high redshift RCS-1 clusters directly shows that BgcR 
over-predicts mass by ~1.4 by z=0.8 compared to z=0.2 (Hicks 2005), though with significant error 
bars. 

The evolution derived in the self-calibration analysis, 0.40^3 80, encompasses these results, due 
to the large uncertainty. In principle we could place priors on the value of 7 with are significantly 
smaller than the formal uncertainty. In practice however, we have left 7 unconstrained; this is 
done because unresolved evolution in the scatter can show up as a large value of 7 (the presence 
of scatter moves the apparent mass limit, and hence evolution in scatter appears as evolution in 
Mum)- We currently have no direct information on the scatter in the mass-richness relation at 
z = 1, and leave 7 unconstrained to account for this. However, note that the fact that our derived 7 
is small implies that the scatter does not evolve strongly. In situ measurements of mass-observables 
in a significant number of clusters at 2; = 1 are needed to progress further; such work is ongoing 
now. Also, in larger surveys, the degeneracy with scatter can also be broken by studying the mass 
function at different redshifts (Lima & Hu 2005). Even modest scatter in the presence of a steep 
mass function results in more up-scatter than down-scatter in mass and so the shape of the mass 
function can be used to calibrate the scatter. 



4.2. Analyses With Mass-Richness Priors 

In order to provide the best possible cosmological constraints from the current cluster sample 

wc have repeated the MCMC analysis using priors on the mass richness relation from both the 
CNOC-1 dynamical results (Yee &: Ellingson 2003) and those of Blindert et al. on an RCS-1 
subsample. The resulting parameter values are summarized in Table 3. The use of external priors 
significantly reduces errors on the cluster parameters, especially those affecting the limiting mass. 
Moreover, external priors make the probability distribution of all the cluster parameters more close 
to Gaussian, but less so for the cosmological parameters. The priors also have the effect of lowering 
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The effect on tfie derived mean value of the cosmological parameters is remarkably small 
however, as might be expected given the degeneracies apparent in Figures 3 and 4. Essentially, 
the degeneracies dominate in the current dataset, and our current observational constraints are 
insufficiently precise to break those degeneracies. We are unwilling currently to constrain the 
scatter and redshift evolution using priors based on available data (though such priors would limit 
the degeneracies currently seen); leaving these parameters free as a check on the consistency of the 
results is a more robust approach. 

5. Conclusions 

The promise of cluster surveys for cosmology has been long held, but in comparison to other 
methods (such as weak lensing, the CMB, or SNe), little realized. The initial general analysis of 
the RCS-1 survey presented here is the first ever attempt to use a large homogeneous sample of 
clusters over a broad mass and redshift range to constrain both fl^ and ag simultaneously. More 
importantly, this cluster sample is both optically selected and characterized; the cluster redshifts 
and mass estimates are derived directly from only two-band imaging to modest depth. It is thus, 
by comparison to many methods, very modest in terms of the required observational resources. 

Using the best available priors on the mass-richness relation, we find nm=0.30t^il and 
O"8=0.70lo'i55 in excellent agreement with the bulk of the literature (and particularly with the 
recent year-three WMAP results). Additionally, in an analysis with no priors on the mass-richness 
relation we find a very similar cosmological result, as well as compelling agreement between our 
constraints on the amplitude, slope, and scatter in the mass-richness relation and equivalent pa- 
rameters derived directly from detailed dynamical studies of a sub-sample of the RCS. Our results 
also compare well with the cluster mass-richness relation from the CNOC-1 cluster sample. The 
success of the analysis presented here, and the consistency of the results both cosmologically and 
with our current understanding of cluster observables, demonstrates that both the RCS-1 cluster 
survey strategy, and the self-calibration methodology are tractable in practice. Perhaps more im- 
portantly it also suggests that clusters, at least in aggregate, are well behaved and amenable to 
being used as cosmological probes to at least z = 1. 

Our self-calibration analysis includes a first effort in treating the scatter in the mass-richness 
relation, using a fixed, mass-independent fractional scatter about the mean mass-richness relation. 
We find the scatter to be significantly degenerate with as ■ These results demonstrate that scatter in 
the mass-observable relation is clearly an important parameter to consider, and dictate that future 
efforts to constrain mass-observable relations (regardless of the observable) probe large enough 
samples with sufficient precision to set significant limits on the scatter, over the entire redshift 
baseline used in such surveys. While this effect may be mitigated using mass observables with 
intrinsically smaller scatter, direct in situ determination of the scatter will remain an important 
component of upcoming large surveys. Detailed follow-up study of a subsample of the RCS-1 cluster 
catalogs with this in mind, as well as the imaging for the more than ten times larger RCS-2 survey. 
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is ongoing; further cosmological results from these data will be presented in future papers. 

Overall, the results shown here hold great promise for much larger optical surveys now under- 
way, ongoing and proposed X-ray surveys, as well as large Sunyaev-Z'eldovich effect-based cluster 
surveys which will begin soon. The analysis presented here is relatively straightforward; additional 
leverage on cosmological constraints in larger surveys can be derived from inclusion of, for example, 
the cluster power-spectrum. Whether or not these additional gains will need to be surrendered to 
a more complex treatment of the mass-observable scatter is not yet clear, but the good agreement 
between our results and current predictions of uncertainties in these large surveys is a positive first 
step towards realizing the promise of galaxy clusters as a precision cosmological probe. 
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Fig. 1. — The input redshift distribution (histogram), along with the best fit cosmological model 
(solid line). Error bars are Poisson only. 
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Fig. 2. — Likelihood functions of the RCS-1 cluster data as function of 6 basic parameters 
used in generating the Markov Chains. There are 3 cosmological parameters (ilmjfS; and h) and 3 
cluster parameters {Mumo,'!, and /sc)- The shaded region shows the marginalized likelihood from 
the Markov chains. The marginalized likelihoods are very non-Gaussian for all cases (as expected 
for cluster counts) except that of h where a Gaussian prior is used. Notice the long tails in fjg and 
fsc- The mean value for each parameter is shown by the solid line. For 0,m and fig we also show 
the 1— and 2 — a regions by the dashed and dotted lines respectively. 
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Fig. 3. — The shaded regions show the joint 1 and 2 a confidence regions for cosmological parameters 
CIm and as and cluster parameters Mumo and 7 which dictates the mass hmit of the survey at 
any redshift. Additional constraints from targeted observations of clusters can break the current 
degeneracies. 
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Fig. 4. — The importance of scatter in RCS-1 constraints. The shaded regions show the joint 1 
and 2 a confidence regions of 0,m, o"8, and Mumg with f^c- The scatter and its dispersion comes 
out of the self-cahbration analysis of the RCS-1 clusters and are in good agreement with targeted 
observations of a subsample of RCS-1 clusters at moderate redshifts. 
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Table 1: Priors in MCMC An alysis. 

Parameter Prior Notes 

0.05 - 0.55 uniform 

as 0.40 - 1.30 uniform 

h 0.72 ± 0.08 Gaussian 

rib 0.046 fixed 

n 0.99 fixed 



-^Bgc 6-14 uniform 

a 0-3 uniform 

7 -4 — 4 uniform 

fsc 0-1 uniform 



Table 2: Derived Parameters from the Self-Calibration Analysis without Mass-Richness Priors. 

Parameter Mean (68% Confidence Range) 

n si+uii 

log{Miimo) 14.61^[]-?^ 
7 OAOtlil 



Table 3: Derived Priramelers from I ho Seir-C'Mlil)rMliou Analysis with Mass-Richness Priors. 





Mean (68% Confidence Range) 


Mean (68% Confidence Range) 


Parameter^ 


(Blindert et al. 2006 Priors) 


(Yee & Ellingson 2003 Priors) 




0.30^-1^ 

0.701°:?^ 


31+"-^^ 
0.68t°;22 


^Bgc 

a 
7 

fsc 


9 61+° *^^ 
0.69l°|° 


1 7^+0.24 
-'-•''-'-0.24 

64+^-^^ 

71+019 
^•'-'^-0.17 



"For the case with no priors, parameter values are as reported in Table 2, with values of Asgc and a as reported in 
the main text (Asgc = lO.SSli:^! and a = IMtoH). 



